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Abstract 

We construct a theory of spin wave excitations in the bilayer manganite 
Lai.2Sri.8Mn207 based on the simplest possible double-exchange model, but 
including leading quantum corrections to the spin wave dispersion and damp- 
ing. Comparison is made with recent inelastic neutron scattering experiments. 
We find that quantum effects account for some part of the measured damping 
of spin waves, but cannot by themselves explain the observed softening of spin 
waves at the zone boundary. Furthermore a doping dependence of the total 
spin wave dispersion and the optical spin wave gap is predicted. 

PACS 75.30.DS, 75.25+z 
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I. INTRODUCTION 



The colossal magnetoresistance (CMR) manganites, of which perhaps the best known is 
Lai_2^Ca^Mn03, have been challenging the theoretical understanding of the way in which 
magnetism and metallic behaviour co-exist for more than fifty years. These materials are 
difficult to describe for precisely the very same reason that they are interesting; namely that 
they exhibit a complex interplay between lattice, charge, orbital and spin degrees of freedom. 
This gives rise to a very rich phase diagram, exhibiting different magnetic, orbital and charge 
orders, and both metallic and insulating behaviour as a function of temperature, pressure, 
applied field and doping iiS. Even within the "simple" low temperature ferromagnetic phase, 
the mechanism for the metal-insulator transition which occurs as ferromagnetic order breaks 
down remains controversial. 

The CMR manganites share a layered perovskite structure with the even more widely 
studied high- temperature (HTc) superconductors; they may be synthesised with one, two, 
three (or many), neighbouring conducting planes. The materials most frequently discussed 
are the three dimensional "infinite layer" compounds, which have equally spaced planes and 
are approximately cubic in symmetry. Here we will construct a theory for ferromagnetism 
in La2-2xSri+2xMn207 and discuss results especially for x=0.4. In this bilayer compound 
planes of magnetic Mn atoms in MnOg octahedra are grouped in well separated pairs. The 
small spin wave dispersion found empirically perpendicular to these planes provides us with a 
justification for considering, as a first approximation, only a single pair of planes i. e. a single 
ferromagnetic bilayer with moments lying in the ab- planeiS. The T-x phase diagram and 
evolution of magnetic structure with doping has been reported iniM and a FM phase persists 
in the range 0.3< x <0.4. For larger doping an intra- bilayer canting of moments appears 
and the charge ordered stoichiometric compound (x=0.5) LaSr2Mn207 finally is an AF 
insulator. Here we will concentrate on predictions for the spin wave dispersion and damping 
of FM bilayer manganites which have been measured by inelastic neutron scatteringBniil. 
Calculation of the spin wave damping requires going beyond the usual semi-classical picture 
used to describe spin wave excitations in the manganites to include quantum effects. In 
Section y we present a minimal model of a bilayer manganite based on Zener's double 
exchange (DE) mechanismlll'El. A fully quantum mechanical large S expansion of this model 
is developed, following a recently introduced operator expansion methodffl. Predictions 
for the dispersion of the optical and acoustic spin wave modes of a double exchange bilayer, 
their doping dependence together with their damping, are made in Section HI, A comparison 
with experimental data for Lai.2Sri 8Mn207 is made in Section 0. This comparison provides 
a test of how well the DE model describes FM in CMR materials when quantum effects are 
included. We conclude in Section |V| with a discussion of the implications of our results for 
the theory of ferromagnetism in CMR manganites. 



II. THE MODEL HAMILTONIAN 

In this section we consider Lai.2Sri.8Mn207, as a concrete example of a bilayer DE system, 
and derive a model Hamiltonian for a single Lai.2Sri.8Mn207 bilayer starting from Zener's 
DE mechanism, in the limit where the strength of the Hund's rule coupling is taken to be 
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infinite. The comparison of tlie predictions of this model with experimental data in Section 
^ therefore provide a test of how well the DE model describes FM in CMR materials. 
The crystal structure of La2_2xSri+2a;Mn207 belongs to the space group I4/mmm with a 
body centred tetragonal conventional unit cell that contains two distorted MnOg octahedra 
as basis whose distortion depends on dopingi'i. The lattice constants are a =3.87 A and 
c= 20.14 A. The intra- bilayer spacing d~a is much smaller than the distance D=6.2 A 
between two adjacent bilayers. Therefore bilayers are well separated, and the spin wave 
spectrum measured by inelastic neutron scattering indeed shows a very small dispersion 
of about 0.4 meV in the direction perpendicular to the planesSS. For this reason we will 
neglect coupling between the planes entirely, and model Lai.2Sri 8Mn207 in terms of a single 
pair of layers. Within a given bilayer, both magnetism and metallic behaviour originate 
in the Mn ci-electrons. Mn t2g (i-orbitals are exactly half filled, and form a spin 3/2 local 
moment because of strong Hund's rule coupling. This local moment couples to itinerant 
electron Cg (i-orbitals through a similar Hund's rule exchange interaction. In the metallic 
phases of the manganites, electrons in Cg orbitals delocalise by hopping between Mn atoms 
through intermediate 02p orbitals — a process named "double exchange" by Zener0. This 
delocalisation of the eg electrons stabilises FM order among the t2g spins, since both are 
strongly coupled by Hund's rule interaction, and the Cg electrons will have the maximum 
kinetic energy if all spins are aligned. 

In the bilayer compounds the MnOg octahedra show a doping dependent pronounced 
Jahn- Teller (JT) distortioni§, therefore Mn^+ site symmetry is no longer cubic and a crys- 
talline electric field (CEF) splitting of e^ (dsz^^r^jd^i-yi) states ensues. The influence of this 
splitting on the stability of magnetic phases was investigated by Okamoto et aP. The e^ 
splitting energy A is generally smaller than the inter-site in-plane hopping t and therefore 
in the FM ground state the orbital state is of uniformly mixed 6.3^2 _^2/dri.2_y2 character. In 
this case orbital degrees of freedom do not appear explicitly in the Hamiltonian but the 
degree of admixture determines the ratio of interlayer (tj,) to intra-layer hopping (t) of the 
effective single band (orbital) Hamiltonian which is then given by 

(ij)Xa ia iXafS 

where c\-^ is the creation operator for an Cg electron on site i of plane A = {1,2} with 
spin a = {T; i}- The components of the operator aa/3 are Pauli matrices, and Si\ is the 
spin operator for the t2g electrons on that site. The on-site exchange Jh parameterises 
Hund's rule coupling, and the sum (ij) runs over nearest neighbours within a plane. Our 
subsequent DE spin wave analysis will lead to t^ 0.175 eV and t_L ~ 0.1 eV. This is much 
smaller than the intra-atomic (Hund's rule) exchange J/^ ~ 2 eV which may be estimated 
from the snlitting of majority and minority spin LDA bands in the stoichiometric (x=0.5) 
compoundBil. 

In addition, there may be super-exchange interactions between spins, both within the 
plane (J) and between the two planes of the bilayer (J±). These can be parameterised by 

Hex = -J'''' ^^-^^-^ - ^^"^ {^^1-^*2 + h. c.} (2) 

(ij)X i 
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FIG. 1. Convention for labelling vertex for interaction between electrons and spin waves in 
the limit Jn/t — > oo. Straight lines correspond to electrons fkrj and wavy lines to spin waves dq^, 
where k and q are momenta in-plane and 77, = 0, vr are the momenta out of plane. 



where Tix = Si\ + l/2J2ai3^i^a(^ai3Cix^ is the total spin operator for both t2g and Cg d- 
electrons on the site iX. Exchange integrals in the manganites can be FM or antiferromag- 
netic (AF) depending on the details of orbital occupancy and electronic structure. 

To evaluate the spectrum, or even to find the ground state of the Hamiltonian Eqn. |l| 
is a formidable task, but if we assume FM order and treat the length of the local moment 
S, and the ratio Jn/t as large parameters, we can derive a controlled expansion of the 
properties of a bilayer ferromagnet. The most direct way of doing this is to work with 
eigenstates of the Hund's rule coupling term, and to quantize small fluctuations of the total 
spin operator Tix using a generalisation of the usual Holstein-Primakoff procedure due to 
Shannon and ChubukoAll. This approach will now be extended to the bilayer system. In 
the limit t/Jn <^ 1, t^/J^ <^ 1 Jh we obtain a model in which bosonic fluctuations of 
the total spin interact with a band of spinless electrons. In this limit it makes sense first 
to diagonalise the Hund's rule coupling term in the Hamiltonian and then to introduce the 
hopping of electrons as a "perturbation" . We do this following the method introduced inil 
by constructing new Fermi operators {/, f^} = 1 and {p,p^} = 1 which create eigenstates of 
the Hund's rule coupling term with eigenvalue —JhS/2 and Jh{S + l)/2, respectively. The 
Hund's rule coupling then reads 



JhS 



/ / - U + o P P+ 



sr s 



(3) 



where the sum over sites has been suppressed. In the physically relevant limit limit Jh ^ 00, 
for less than half filling, we can remove p Fermions from the problem entirely, and rewrite the 
kinetic energy term in Eqn. |l| entirely in terms of a band of spinless (/) electrons interacting 
with fluctuations of the total spin parameterised by the Bose operators [a, = 1. 

To accomplish this transcription of the Hamiltonian it is sufficient to know a few of the 
leading terms of the inverse transformation between C| and Cj electron operators, and the 
new / operators creating eigenstates of the Hund's rule coupling term 

c,^/(l-§).... (4) 
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FIG. 2. Leading self energy corrections for spin waves due to interaction with electrons in 
limit Jh/t oo. Diagrams Ia~f) show contributions for acoustic (flgo) ^-nd diagrams Ila-f) optical 
(ogTr spin wave modes. 
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To prove this result, and to derive the full transformation between "laboratory frame" C| and 
Cj electron operators and the "local frame" / and p operators, together with the appropriate 
algebra for the spin boson a is an involved task. We will not discuss the transformation in 
detail here (seeEj), but note that all the necessary canonical commutation and anticommu- 
tation relations, e. g. [/, a] = 0, etc. , are obeyed. 

Up to a constant the transformed DE Hamiltonian reads 



n = no + V2 + o{i/s^) 

where the kinetic energy term for / electrons is 

no = -tJ2 flfjx -t±Yl {f^^^ + h.c.] 



(6) 



(7) 



At this level the spin excitations a are dispersion-less. Spin wave dispersion first enters 
into the problem through interaction at 0{1/S) through the interaction term 



+ h. c. \ (8) 



where we have neglected a further four boson vertex at 0{1/S'^) which is irrelevant at zero 
temperature. 
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By Fourier transformation we obtain the following Hamiltonian which describes a band 
of spinless electrons interacting with (initially dispersion-less) bosonic spin-wave excitations. 

n = no + V2 + o{i/s^) 

no = Yl - t±) fUko + {ek + t±) flfk. (9) 

k 

fci---fc4 ■rjrj'i^i/'={0,TT} 

where we consider symmetric and antisymmetric combinations of electron operators (binding 
and antibinding bands), and of spin operators (acoustic and optical spin waves), for the two 
planes. For the simple nearest neighbour tight-binding model Eqn. |l] the in-plane electronic 
dispersion is given by = - zt^( cos kx + cos fcy) in units where the distance between Mn 
atoms a = 1. The scale of interaction between electrons and spin waves V2 is determined 
entirely by electronic energies, but is one order in S down on the kinetic energy term Hq. 
There are a total of eight physically distinct vertices (decay channels) for interaction between 
electrons and spin excitations. The convention for labelling these vertices is shown in Fig. |I| 
and their algebraic expressions are given in Eqn. and [A2| . The spin wave dispersion is 
now determined by the leading order self energy up to 1/S^ shown in Fig. We first discuss 
the results within the usual semiclassical (1/S) approximation. 



III. THEORETICAL PREDICTIONS 



In a cubic system, at a semi-classical level of approximation, Zener's DE mechanism leads 
to a FM effective nearest neighbour Heisenberg exchange interaction between neighbouring 
Mn spins, with a spin wave dispersion 

o;, = zJ^^5[l-7,] (10) 
where the size of the effective exchange interaction is set by electron energiesEl 

k 

Here ■jg = ^{cosqx + cos qy + cosg^) is the structure factor for a 3D cubic lattice and n{k) 
is the occupation of the electronic state with momentum k, and J^^ is proportional to the 
expectation value of the kinetic energy operator per Mn-Mn bond, relative to the center of 
the band. Spin waves are exact eigenstates of a Heisenberg FM, and therefore undamped. 
This simple mapping between DE and Heisenberg models breaks down, however, when 
quantum effects are included@. 

The situation in a bilayer system is complicated in that there are both optical and 
acoustic branches of spin wave excitations, but the mapping onto an effective Heisenberg 
model is once again possible at a semi-classical level. Evaluating the effect of interaction 
between electrons and spin waves described by Eqn. || to 0{1/S), and now including the 
effect of super-exchange terms, we obtain a spectrum : 
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FIG. 3. Inset shows total DOS Nt(e)= No(e)+ N^(e) and total electron number n for the 2D 
bands eo{k) and e,r(^) Wt=W+t^ is half of the overall bandwidth. Main figure shows variation of 
and Jj?^ = with hole doping under the assumption that t_L/t=0.57 (fixed for x=0.4) is 
independent of x. The physical FM regime is restricted to 0.3< x <0.4. 



= z{j'''' + J^^)5[l - 7j + 2( ^ + Jf (12) 



1 

where in 2D 7^ = ^(cosg^. + cosgj^), of*^ is the dispersion of the acoustic and cuj the dispersion 
of the optical spin wave branch. The size of the DE in-plane contribution to the effective 
exchange integral is once again set by the expectation value of the kinetic energy on a single 
bond, and the DE between the two planes is determined by the occupation difference of 
binding and antibinding bands: 

k 

'^^''=2^^S[''°^^^~'''^^^^] ^^^^ 

k 

Here we used the occupation numbers n^ik^ ={flofko) and nT^(k) ={f\T^fk-K) of the binding 
and antibinding electron bands eQ{k)= -t_L+e(A;), e.„{k)= t_L+e(fc) respectively. Our result at 
this order agrees perfectly with earlier calculations of the spin wave spectrum in a bilayeiSS. 
The effective exchange constants in Eqn. |13| can be evaluated as function of the doping x 
which gives the number of holes per Mn-site or the total number of e^ electrons per Mn site 
n= no+n7r= 1-x that occupy the 0,7r- bands. By using the DOS functions No,7r(e)=N(e±t_|_), 
the electron number uq^tt and the average band energy eo.Tr of the 2D binding and antibinding 
bands respectively may be expressed as 
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FIG. 4. Doping dependence of DE exchange constants. The strong reduction of jf ^ =J^^ 
with increasing x is due to the strong reduction of tc =t± due to__the JT effect of MnO units 



as described by Eqn. 17 with Q±r]±=-2. Experimental data fromllSl (circles) Ji (squares) ani 
(triangles). In the latter case a somewhat smaller SJ^^=8.6 meV was obtained by fitting the spin 
wave stiffness constant in contrast to the value of 10 meV obtained by fitting to the whole spin 
wave bandwidth will?''^=zSJ^'^. The inset shows the dependence of DE exchange anisotropy on 
the hopping anisotropy for a fixed doping x=0.4. 



N{e)de 

eo,. = r^'" N{e)ede (14) 
J-w 

Here W=zt and 2W is the band width of of each of the 2D bands eo,7r(^) and ei;' is the Fermi 
leveL Furthermore K(,^) is the complete elliptic integral of the first kind. The total DOS 
Nt=No+N.n- and the total number of electrons n as a function of the Fermi level is shown 
in the inset of Fig. The spikes in the DOS are logarithmic singularities of each of the 
2D bands at its band center (itj,). We then obtain for the effective DE exchange constants 
after Eqn. ^ 

JS" = ("0 - "J (15) 

The DE anisotropy ratio Jf^/J^^ is is equal to the ratio cuj/wi^^' of optical spin wave 
gap ci;q=2SJ;^^ to the acoustical (or optical) spin wave band width wili°^= 2zSJ^^ along 
[110] direction and it is given by 
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J^^ ^t^ eo + e^ ^ ^ 

For t± <^t the occupation number difference increases hnearly in t^, and using Eqn. ^ 
we find J^^/J^^ ~ {t±/t)'^, as shown in the inset of Fig. ^. Numerical values from Eqs. 



151 are presented in Fig. ^ as function of the doping x=l-n. It shows the variation for the 
model DE for fixed t^ in the whole range 0< x < 1 although it must be kept in mind that 
the physical region for the FM phase of La2_2a;Sri+2xMn207 is much smaller, according tou 
it exists for 0.3< x < 0.4. From a comparison of the experimental values of the optical 
spin wave gap and the spin wave band width at x=0.4 with Eqn. ^ and with the insert in 
Fig. ^ we can obtain estimates for the underlying microscopic model parameters within the 
classical approximation, namely tx/t~ 0.57 corresponding to the experimental J^^/J^^ ~ 
0.30 at low temperature and t~ 0.175 eV (t^= 0.1 eV) as obtained from the experimental 
value SJ'^'^= 10 meV (from w[j^°'= zSJ^^= 40 meV0) by using Eqn. |15[ According to Fig. 
^ J^^(x) and J^^(x) should not change dramatically with the hole doping in the FM regime 
0.3 < X < 0.4, namely at most 6% and 15% respectively. However this refers to the artificial 
situation where t_L does not depend on the doping. From the experimental investigation of 
optical spin wave gap and dispersion for various dopings x=0.30, x=0.35 and x=0.40li3 it is 
known that indeed J^^ shows no change in this region, however strongly increases 

by a factor of four when the doping is reduced from x=0.4 to x=0.3. 

The origin of this pronounced doping dependence of interlayer DE is connected to the 
large Jahn-Teller(JT) distortion observedEl in the bilayer manganites. This distortion is 
defined as AjT=apical Mn-0 bond length/equatorial bond length. Decreasing the doping 
leads to an increase of Ajt- The driving mechanism for this JT distortion is an increasing 
admixture of d^z^-r^ states into the conduction band states which naturally leads to an 
increase of t^ with reduced doping, which in turn strongly increases the interlayer Jx^(x) 
as shown in the inset of Fig. ^. For the limited FM doping range one may describe this de- 
pendence by introducing dimensionless Griineisenparameters T]±=-{d\n D/dlnx) and f2x=- 
{dlnt±/d\nD) where D= distance between the layers of a single bilayer. They describe the 
doping dependence of the JT- distortion and the distortion dependence of the interlayer- 
hopping respectively. The JT effect on the intra-layer hopping t is neglected since no doping 
dependence of J^^ is observed. Assuming that ri± and Q± are constants in the range of x 
considered, this amounts to a doping dependence of t^ given by 

t±(x)=t°(-)^-''- (17) 

Xo 

where e.g. xo=0.4 and t'[=tx(xo). According to the physical origin of the JT distortion 
mentioned above one has to expect that ^l±'r]± < 0. Using the above relation in Eqs. ^ 
with fi_L?7_L= -2 one obtains the doping dependence of the exchange constants shown in Fig. 
^ together with the experimental values for various dopings. Using ri± ~ 0.037 from the JT- 
distortions given ini one then obtains from the above relation fi^ ~ -54. This large negative 
Griineisenparameter characterises a strong dependence of the effective interlayer hopping 
t^ on the layer spacing D within a bilayer. The JT- distortion increases with temperature 
for constant doping implying an increase of t^ and hence J^^. The DE therefore becomes 
more isotropic at higher temperature. This has indeed been observed for x=0.4 in diffuse 



9 



Spinwave 
Energy [eV] 0.1 




FIG. 5. Spectrum for optical and acoustic spin wave modes throughout the Brillouin zone, 
calculated for doping x = 0.4, t = 0.175ey, t± = O.leV. Upper pair of lines — semi-classical spin 
wave dispersion at 0{1/S); lower pair of lines — spin wave spectrum including quantum effects at 
0(1/52). 



neutron scattering where J^^/J^^ ~ 0.5 has been found at room temperatureE3 compared 
to ^ 0.3 from the low temperature spin wave experiments discussed here@. 

The results for the dispersion of acoustic and optic spin wave modes for throughout the 
Brillouin zone, at a semi-classical level and for the parameters given above are shown in 
Fig. ^ At zero temperature, at a semi-classical level, spin waves are undamped, i. e. states 
with a single spin wave excitation are good eigenstates of the Hamiltonian, with no allowed 
decay processes. The dispersion of spin waves is generated by their elastic scattering by the 
average density of electrons. At 0{1/S'^) spin waves can decay through interaction with 
electrons into lower energy spin excitations dressed with particle hole pairs — an inelastic 
process. This leads to a damping of spin waves, and a corresponding shift in spin wave 
dispersion to lower energy. We can evaluate both effects starting from the Hamiltonian 
Eqn. |[ We consider 

^,^ = o;,^ + i?e{S^(g,0)} 

7,^ = -/m{S'^(g,^P} (18) 

where a;^ is the net dispersion and 7^ the damping of the spin wave excitation, and S'^(g, Q) 
is the momentum and frequency dependent self energy correction due to interaction of spin 
waves with electrons at (^(l/S'^). The various contributions to the spin wave self energy are 
shown in Fig. |^ at this order and given in Appendix |^. The new physical process involved 
is the inelastic scattering of spin waves from fluctuations of charge density. Results for spin 
wave dispersion, including leading quantum corrections are shown in Fig. ^, and values for 
the damping of spin waves in Fig. 
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FIG. 6. Damping of acoustic and optical spin wave modes, throughout the BriUouin Zone, 
calculated for doping x = 0.4, t = 0.175ey, t± = Q.leV . 

The immediate conclusion which we can draw from these calculations is that quantum 
effects on spin waves in a DE bilayer are very large. The downward renormalisation of spin 
wave dispersion at Oil/S"^) is a sizable fraction (about 30%) of the spin wave dispersion 
at 0{1/S). Similarly, the damping of spin waves is quite pronounced, being of the scale 
5-10% of the spin wave dispersion, rising to a maximum value ~ 6 meV at the zone corner. 
Because of the large renormalisation of the spin wave spectrum at 0{1/S'^) it would be 
necessary to reparameterise our model to fit experimental data with the leading quantum 
effects included, by increasing the sizes of the hopping integrals t and t_L, and including 
super-exchange interactions J^^ and J_f^. Any increase in the electron bandwidth would 
give a proportionate increase in the damping of spin waves. 

Examining the quantum corrections in more detail, we find that the spin wave dispersion 
has been modified so as to give a relative softening of spin wave modes near the zone center. 
This can be understood loosely in terms of the dynamical generation of an effective non 
nearest-neighbour couplings between spins by processes Oil/S"^). It is also interesting to 
note that the gap between acoustic and optical modes is now momentum dependent. 

While these effects are of themselves interesting, they do not offer any unambiguous 
signatures of quantum effects in the magnetism of La2_2xSri+2a;Mn207, as one could achieve 
similar modifications of the spin wave dispersion simply by postulating additional exchange 
couplings between spins on an ad hoc basis. It is the damping of spin waves at zero tempera- 
ture which sharply distinguishes a DE system from any conceivable Heisenberg ferromagnet. 
In a the Heisenberg ferromagnet, spin waves are undamped at zero temperature, and damp- 
ing only becomes appreciable for temperatures large compared with the spin stiffness D. 

The damping which we predict for the DE model Eqn. |I| at zero temperature is large 
and highly momentum dependent. The zone centre acoustic mode must remain undamped 
(in the absence of any magnetic anisotropy) as it corresponds to the rotation of the total 
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magnetization of the is the Goldstone mode of the system. Accordingly, in the zone centre, 
we find that the damping of the acoustic mode vanishes as 

Tf = a^^q' (19) 

The optical spin waves are not Goldstone modes, however, and have a finite dispersion (and 
damping) in the zone centre. We find that the latter behaves as 

= + a^^g^ (20) 

The lower power law in q here reflects the way in which the vertex for spin wave scattering 
is cut off by the interlayer hopping t±. 

Away from the zone centre the spin wave damping exhibits stationary points at the 
symmetry points of the brillouin zone — a maximum for both acoustic and optic modes 
at X, and a minimum for both at M. It is interesting to note that the higher energy 
optic modes are not always more strongly damped than the acoustic modes, and that the 
maximum damping does not occur for the highest spin wave energies, as one might expect. 
Infact the momentum dependence of the damping of spin waves in DE systems varies strongly 
with doping, being constrained by both the geometry of the Fermi surface and the complex 
momentum dependence of the spin wave scattering vertex. 



IV. COMPARISON WITH EXPERIMENT 

The spin wave dispersions in bilayer manganites have been investigated by several 
groupsocJ, and a consensus was reached that the data could not be explained using a 
nearest neighbour Heisenberg model dispersion of the form Eqn. Departures from 

Heisenberg model behaviour also have been observed many cubic manganite systems, for 
example Pro.eaSro.sTMnOacSl. Typically, what has been seen in both cubic and bilayer sys- 
tems is a softening and broadening of the zone boundary spin waves. The total spin wave 
bandwidth measured to the zone boundary is much less than would be predicted on the basis 
of the spin stiffness D measured in the zone center, and the zone boundary spin wave modes 
are extremely broad in comparison with their energy. The theory of bilayer manganites pre- 
sented here shows that the double exchange model can exhibit both of these effects, when 
quantum corrections are included. However, the minimal model Eqn. |l], as parameterized 
above, is not sufficient to obtain a quantitative description of the experimental results. 

With regard to the spin wave dispersion, the inclusion of the quantum corrections shown 
in Fig 1^ cannot explain the measured departures from Eqn. |12|. The dispersion measured by 
inelastic neutron scattering on Lai.2Sri 8Mn2O70, for the acoustic mode in the zone center, 
has the form 

u;? = A + (21) 

with A ^ D^, as would be expected for a FM with small magnetic anisotropy. However at 
larger momentum transfer the measured dispersion lies below the curve 

u'^ = zD'[l - (22) 
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(away from the zone center we can safely neglect A). The total acoustic spin wave band- 
width, as defined by the spin wave energy at the zone corner, is about 15% less than 2zD. 
If we compare with a suitably parameterized quantum theory for the DE model, we find 
that the predicted dispersion lies above the curve Eqn. and so the observed softening 
effect is absent — infact quantum corrections have the wrong "sign". This failure of the 
minimal model Eqn. |I] is neither very surprising nor very disappointing, given that we have 
attempted to fit the spin wave dispersion of a complex system with spin charge and lattice 
degrees of freedom throughout the Brillouin Zone, using only two adjustable parameters. 
However it is important to ask which of the many simplifications made is to blame for this 
disagreement with experiment ? 

A better fit could probably be obtained at a semi-classical level, by substituting a more 
realistic dispersion for the underlying electrons into the one loop diagrams used to calculate 
the 0{1/S) spin wave self energy. In tight binding language, each hopping integral tij has 
a corresponding DE coupling J^*^ associated with it. The inclusion of tij beyond nearest 
neighbours to obtain a more realistic electronic bandstructure therefore also modifies the 
form of dispersion of the classically equivalent effective Heisenberg model. Attempts to 
calculate spin wave dispersion directly from electronic structure suggest that this effect is 
important, and leads to a softening of zone boundary modes, at least in cubic systems@. 

At a quantum mechanical level, since interactions between spin waves are mediated by 
density fluctuations of the electron gas, it would be more realistic to use a screened form of 
the charge susceptibility in which long range interactions were suppressed. We anticipate 
that this would also tend to enhance the softening of zone boundary modes. The inclusion 
of leading Quantum corrections in 0(t/ Jh), likewise leads to a softening of zone boundary 
spin wavesil'0. 

Each of these improvements to the model would involve the introduction of new param- 
eters, which would need to be checked against electronic structure and other experiments. 
Since the stated aim of this paper is to explore the minimal model Eqn. |I|, we will not discuss 
such refinements further here. A more interesting possibility to explain the difference be- 
tween experiment and theory would be that spin waves are coupled to orbital and/or lattice 
modes. We return to this below. 

The present theory is more successful in explaining the damping at least in a qualitative 
way. Fig. |^ shows the experimentally observed widths of the energy scans as a function of 
momentum transfer q from zone center to the zone boundary for the acoustic spin excitations 
of bilayer manganite Lai.2Sri 8Mn207, along with that obtained from the present theory 
(continuous line). Data measured at different Neutron sources have been plotted together. 

The damping which we calculate has a similar momentum dependence to that observed, 
but is smaller by approximately a factor of four. Some part of the difference in absolute value 
between theory and experiment can be explained by the fact that the parameters t and t± 
were chosen so as to correctly reproduce the spin wave bandwidth at a semi-classical level, 
and somewhat larger values must be used to fit the measured dispersion once quantum effects 
are included, leading to some systematic variation in the values quoted and their associated 
errors. 

Our conclusion is that the minimal DE model fails to explain the softening of zone 
boundary spin waves in Lai.2Sri.8Mn207, but can explain about 30-40% of their width. Of 
the refinements to the model discussed above, only the use of a screened charge susceptibility 
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FIG. 7. Damping of acoustic spin wave modes on the line TX. Points with error bars — exper- 
imental data for Lai.2Sri.8Mn207 taken fromi; lower solid line — theoretical predicted damping 
for minimal model with x = 0.4, t = 0.175eV , t = O.leV. The dashed line is a guide to the eye. 



would affect the calculated damping of spin waves, but we do not anticipate that this would 
lead to a marked increase in their width. We therefore conclude that spin waves in the 
FM phase of this bilayer manganite are coupled dynamically to another mode, probably of 
orbital and/or lattice excitations. Such a couplings have been proposed in the context of the 
cubic manganites — for exarnple to optical phonon^^l or through Jahn Teller active lattice 
modes to orbital fluctuationsEB. 



V. CONCLUSIONS 

We have constructed the simplest possible model for ferromagnetism in Lai.2Sri.8Mn207, 
based on Zener's double exchange mechanism within a one orbital picture for a single bilayer. 
This model has two adjustable parameters, the intra- and inter-plane hopping integrals t 
and t±. At a semi-classical level it is equivalent to a Heisenberg model with intra- and 
interplane exchange integral J^^ and J^^. The doping dependance of these parameters 
was discussed, and the predictions of the effective Heisenberg model compared with the 
results of inelastic Neutron scattering experiments. 

As the experiments show departures from simple Heisenberg model behaviour in both the 
form of dispersion and the scale of damping of the spin waves at low temperatures, we also 
calculated the leading quantum corrections to spin wave self energies. These arise because of 
the scattering of spin waves from density fluctuations of the electron gas which are neglected 
in the semi-classical approximation. We find that the minimal model considered cannot 
explain the softening of zone boundary spin wave modes, and somewhat underestimates the 
damping of spin waves, even when quantum corrections are included. This suggests that 
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spin waves are strongly coupled to another low energy mode, presumably related to lattice 
fluctuations, either by a direct coupling to phonons or an indirect "orbital fluctuation" effect. 
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APPENDIX A: SPINWAVE ELECTRON INTERACTION VERTICES AND SELF 

ENERGY CORRECTIONS AT 0[l/S'^) 



First we give the interaction vertices ^^^l^i in Fig. |1| and Eqn. |^. There are eight possible 
channels for electron-spin wave interaction; these are labeled according to the convention in 
Fig. The coefficients of these vertices are given by 
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Where the vertex depends on in-plane momenta only through in-plane electronic dispersion 

(A2) 
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where = —zt^{cos kx + cos ky). The fundamental energy scales in the DEFM are set by 
the kinetic energy of the itinerant electrons, and so it is natural that the electron spin wave 
scattering vertices are proportional to t/t_i. 

Knowledge of the Hamiltonian Eqn. ^ is sufficient to develop a zero temperature dia- 
grammatic perturbation theory in 1/5* for the spin wave dispersion of the DE bilayer, up 
to (9(1/5'^), and to calculate the leading contributions to spin wave damping. The relevant 
processes are shown in Fig. |^. At 0{1/S) only the single loop diagrams a) and b) contribute. 
These evaluate to give the Heisenberg-model like result Eqn. |12| for the semi-classical spin 
wave dispersion. 

The one loop diagrams also contribute a constant term and a further renormalisation 
of the classical dispersion at (9(1/5'^), but all new quantum effects arise from the new 
processes contributing to spin wave self energy at the "watermelon" dia grams 

shown in Figure H c)-f). The self energy corrections for acoustic modes evaluate to give : 
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where we have written out electron energies expUcitly and suppressed terms of 0{1/S^) in 
the vertex. The corresponding processes for optical spin waves yield : 
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To C( 1/5*2), we can neglect the frequency dependence of the denominator in these expres- 
sions and evaluate the leading quantum corrections to the dispersion of optical and acoustic 
spin wave branches numerically by Monte Carlo integration. 

If we restore the frequency dependence of the self energy terms, we can also calculate 
the imaginary part of each. We can use this to estimate the spin wave damping on the mass 
shell, by setting the external frequency equal to the semi-classical spin wave dispersion at 
that wave number, i. e. setting Q = <^°''^, and eliminating all terms in the numerator of order 
spin wave frequencies. The contribution to damping from diagram lid) is, for example : 
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